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Foreword 


A  knowledge  of  ice  thickness,  ice  motion,  and  ice  concentration  (ice  edge) 
is  essential  for  effective,  safe  ship  operations  in  polar  regions.  The  time 
evolution  of  these  parameters  can  impact  the  performance  of  weapons  systems, 
acoustic  surveillance  capabilities,  search  and  rescue  planning,  and  other  aspects 
of  naval  operations.  The  Naval  Ocean  Research  and  Development  Activity 
has  developed  the  first' high-resolution  sea-ice  forecasting  model  for  the 
Barents  Sea  based  on  the  Hibler  ice  model.  The  model  will  run  operationally 
at  the  Fleet  Numerical  Oceanography  Center  using  ice  thickness  fields  from 
the  Polar  Ice  Prediction  System  as  boundary  conditions.  Testing  has  shown 
that  the  high-resolution  models  are  capable  of  better  forecasting  of  the  ice 
edge  location  and  the  growth  (decay)  and  movement  of  ice  near  land 
boundaries.  The  operational  version  of  the  Barents  Sea  model  will  produce 
the  same  14  sea-ice  products  for  the  Barents  Sea,  on  a  daily  basis,  that  are 
presently  produced  by  the  Polar  Ice  Prediction  System  for  the  Arctic. 


W.  B.  Moseley 
Technical  Director 


J.  B.  Tupaz,  Capt,  in,  USN 
Commanding  Officer 


Executive  Summary 


The  hydrodynamic/thermodynamic  Arctic  sea-ice  mode)  designed  by 
W.  D.  Hibler  of  the  Cold  Regions  Research  and  Engineering  L.aboratory 
(CRREL)  has  been  adapted  to  the  Barents  Sea.  This  model  is  driven  by 
atmospheric  forcing  from  the  Naval  Operational  Global  Atmospheric 

•  Prediction  System  (NOGAPS)  and  oceanic  forcing  from  the  Hibler  and  Rrvan 

coupled  ice-ocean  model.  This  high-resolution  model  (25  km),  which  covers 
the  entire  Barents  Sea  and  the  western  half  of  the  Kara  Sea,  uses  a  6-hour 
time  step.  Development  of  this  model  required  the  design  of  new  ice 
inflow/outflow  boundary  conditions,  which  use  the  ice  thickness  fields  from 
the  Polar  Ico  Prediction  System  (PIPS)  when  inflow  is  indicated.  Model  results 

•  show  good  agreement  with  such  data  as  the  Naval  Polar  Oceanography  Center’s 

(NPOC)  analysis  of  ice  concentration  and  concentrations  derived  from  passive 
microwave  data.  The  model  has  a  tendency,  however,  to  melt  ice  too  quickly 
in  summer  and  to  grow  it  back  too  slowly  in  the  fall.  Planned  improvements 
in  the  atmospheric  and  oceanic  forcing  should  correct  this  problem.  The  high 
resolution  of  the  Barents  Sea  model  enables  it  to  predict  the  ice  edge,  ice  growth 

•  and  decay,  and  the  movement  of  kx  near  land  boundaries  with  greater  accuracy 

than  does  the  PIPS  model.  The  Regional  Polar  Ice  Prediction  System  for  the 
Barents  Sea  (RPIPS-B)  is  the  forecast  system  designed  to  run  at  the  Fleet 
Numerical  Oceanography  Center  (FNOC)  based  on  the  Barents  Sea  ice  model. 
RPIPS-B  is  updated  weekly  by  the  NPOC  analysis  of  ice  concentration.  The 
forecast  system,  presently  in  its  “operational  checkout”  phase,  is  being  made 

®  ready  for  a  winter-spring  operational  test. 
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The  Regional  Polar  Ice  Prediction  System — Barents  Sea  (RPIPS-  B): 
A  Technical  Description 


I.  Introduction 


The  capability  to  forecast  sea-ice  characteristics  in 
the  world’s  polar  oceans  has  become  a  more  viable  task 
in  recent  years.  The  existence  of  increasingly  more 
powerful  computers,  larger  quantities  of  remotely 
sensed  data,  and  the  development  of  improved  sea-ice 
models  has  contributed  to  the  design  of  accurate 
sea-ice  forecasting  systems. 

In  September  1987,  the  Navy’s  Polar  Ice  Prediction 
System  (PIPS)  became  an  operational  model  at  the 
Fleet  Numerical  Oceanography  Center  (FNOC).  The 
model  covers  the  central  Arctic,  the  Barents  Sea,  and 
northern  half  of  the  Greenland  and  Norwegian  Seas. 
PIPS  is  based  on  the  Hibler  ice  model  (Hibler,  1979; 
1980).  The  ice  model  is  driven  by  ocean  currents  and 
ocean  heat  fluxes  derived  from  the  Hibler  and  Bryan 
ice-ocean  model  (Hibler  and  Bryan,  1984;  1987). 
Atmospheric  forcing  for  the  model  comes  from  the 
Navy’s  operational  atmospheric  forecast  model,  the 
Naval  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  (Rosmond,  1981).  PIPS  makes  a 
120-hour  forecast  each  day  of  such  sea-ice 
characteristics  as  ice  drift,  ice  thickness,  and  ice  con¬ 
centration  (ice  edge).  The  length  of  this  forecast  is 
based  on  the  length  of  the  NOGAPS  forecast.  The 
model  is  updated  once  a  week  by  a  digitized  analysis 
field  of  ice  concentration  derived  by  the  Nasal  Polar 
Oceanography  Center  (NPOC).  On  days  when  the 
update  is  not  available,  the  model  is  restarted  from  its 
own  24-hour  forecast  fields  (see  Preller,  1985;  Preller 
and  Posey,  1989). 

PIPS  was  originally  designed  to  cover  a  large  portion 
of  the  northern,  ice-covered  polar  oceans.  The 
resolution  of  PIPS  was  chosen  to  minimize  the  amount 
of  computing  time  used  in  an  operational  forecast. 
Once  PIPS  became  operational,  higher  resolution 
models  designed  to  cover  only  one  particular  region 
of  the  Arctic  would  be  developed  and  run 
operationally.  The  advantages  of  these  high-resolution 
regional  models  are  a  more  accurate  prediction  of  the 
ice-edge  location;  greater  accuracy  in  defining  land 
boundaries,  including  island  boundaries;  and  better 
resolution  of  straits. 

The  first  such  regional  model,  designed  by  the  Naval 
Ocean  Research  and  Development  Activity  (NORDA), 


is  a  model  for  the  Barents  Sea.  The  Barents  Sea 
model  is  similar  to  PIPS  in  design.  The  model  uses  the 
Hibler  ice  model  as  its  basis  and  is  driven  by  NOGAPS 
atmospheric  forcing  and  the  Hibler  and  Bryan 
ice-ocean  model’s  oceanic  forcing.  The  main  difference 
between  the  two  models  is  in  the  boundary  conditions. 
PIPS  uses  solid-wall,  closed  boundaries  everywhere 
except  at  the  southern  boundary  hetween  Greenland 
and  Iceland  and  Iceland  and  Norway.  T  his  boundary 
can  be  simply  treated  as  an  ice  outflow  boundary.  The 
Barents  Sea  model,  however,  required  the  development 
of  new  “ice”  inflow/outflow  boundary  conditions. 

This  report  provides  a  technical  description  of  the 
Barents  Sea  model  and  presents  results  from  both  its 
development  and  operational  checkout  phases.  The 
operational  system  for  the  Barents  Sea  will  be  referred 
to  as  the  “regional  PIPS  for  the  Barents  Sea,”  or 
RPIPS-B. 

II.  Model  Description 

The  regional  Barents  Sea  model  uses  the  Hibler 
dynamic/thermodynamic  sea-ice  model  as  its  basis.  The 
ice  model  is  defined  by  five  major  components: 
momentum  balance,  ice  rheology,  ice  thickness 
distribution,  ice  strength,  and  air  ice/ocean  heat 
balance. 

Although  this  model  description  is  found  in  Preller 
and  Posey  (1989),  it  will  be  included  in  this  report  for 
the  sake  of  completeness. 

The  momentum  balance  used  to  determine  ice  drift 
is  given  by 

Du  »—*■—*■—*  ■ 

m  —  =  mfkxu  +  Ta  +  r„  -  mg  grad  H  +  F, 

where  m  is  the  ice  mass  per  unit  area,  u  is  the  ice 
velocity, /is  the  Coriolis  parameter,  tu  and  r„  are  the 
air  and  water  stresses,  g  is  the  acceleration  of  gravity, 
H  is  the  sea  surface  dynamic  height  and  F \s  the  force 
due  to  variation  in  the  internal  ice  stress.  Ice  is 
considered  to  move  in  a  two-dimensional  field  with 
forcing  applied  through  simple  planetary  boundary- 
layer  formulations. 

The  air  and  water  stresses  are  defined  using  constant 
turning  angles 
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The  equations  for  thickness  and  compactness  arc 


-  Pu  Cd  L\  |  (i.\  cos  <)>  +  k  x  Hk  sinty) 

T„  =  p„  CK  |  Un.  -  u  \  [ (UK  -  u )  cosQ 
+  k  x  (Uw  -  u )  sind  j  , 

where  u  is  the  ice  drift  velocity,  U  is  the  gcostrophic 
wind,  liH  is  the  gcostrophic  ocean  current,  Cu  and  C„ 
are  the  air  and  water  drag  coefficients,  ou  and  pn  are 
the  air  and  water  densities  and  <j>  and  9  are  the  air  and 
water  turning  angles.  For  a  more  detailed  discussion 
of  model  dynamics  and  the  spatial  finite  differencing 
code,  see  Hibler  (1979). 

The  ice  rheology,  a  viscous-plastic  constitutive  law, 
relates  the  ice  stress  to  ice  deformation  and  ice  strength 
in  the  following  manner: 

=  2r]  (t,r  P)  +  lL(t,r  P)  -  n  {£„,  PH 
(5„  -  Pd,/ 2, 


where  a  is  the  two-dimensional  stress  tensor,  tu  is  the 
strain  tensor,  P/2  is  a  pressure  term,  and  £  and  r\  are 
nonlinear  bulk  and  shear  viscosities.  Ice  flows 
plastically  for  normal  strain  rates  and  deforms  in  a 
linear  viscous  manner  for  small  strain  rates. 

The  ice  thickness  distribution  takes  into  account  the 
ice  thickness  evolution  as  a  result  of  dynamic  and 
thermodynamic  effects.  The  regional  Barents  Sea 
model  originally  used  a  two-level  approach  (Hibler, 
1979).  This  approach  breaks  ice  into  two  categories, 
thick  and  thin,  with  the  division  between  the  two  being 
0.5  m.  The  compactness,  .9,  is  defined  as  the  area 
within  a  grid  cell  covered  by  thick  ice,  while  (I  -  A) 
is  the  area  covered  by  thin  ice.  This  treatment  resulted 
in  an  average  ice  thickness  over  the  Arctic  that  was 
too  thin  when  compared  to  observations  (Preller  et  al., 
1986).  I  o  correct  for  this  bias  and  to  include  the  strong 
dependence  of  ice  growth  rates  on  thickness,  a  seven- 
level  ice  thickness  calculation  used  by  Walsh  et  al. 
(1985)  was  added  to  the  PIPS  model  in  March  1988. 
This  method  divides  the  “thick”  ice  into  seven 
categories  and  allows  ice  to  grow/decay  in  each 
category.  The  seven  levels  are  equally  spaced  between 
0  and  twice  heff,  where  heff  is  the  effective  ice  thickness 
(Hibler,  1979)  or  mean  ice  thickness  over  the  entire  grid 
cell.  For  periods  of  ice  growth,  snow  cover  is  also 
divided  into  a  seven-level  linear  distribution  of  snow 
depths  equally  spaced  between  0  and  2  times  the  grid 
cell  mean.  W  hen  melting  occurs,  snow  is  assumed  to 
be  uniformly  distributed  over  the  ice  covered  portion 
of  the  grid  cell.  Snowfall  rates  are  based  on  monthly 
mean  climatological  values  (Mavkut  and  Untersteiner. 
1969;  Parkinson  and  W  ashington,  1979).  W  hen  tested 
on  the  PIPS  model,  this  improved  treatment  of  thick 
ice  resulted  in  an  average  increase  of  ice  thickness  of 
50  cm  over  the  PIPS  domain. 
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where  Sh  and  S,  are  thermodynamic  terms  defined  by 
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with  f(h)  as  the  growth  rate  of  ice  of  thickness  h  (heff) 
and  ha  a  fixed  demarcation  between  thick  and  thin 
ice.  In  all  model  simulations,  h0  =  0.5  m.  In  the 
seven-level  ice  thickness  calculation,  heff  is  the  seven- 
level  sum  of  ice  thickness,  including  the  calculated  snow' 
and  ice  thickness  changes.  The  term  Sh  is  the  net 
growth  or  melt  of  ice.  S.,  is  the  change  in 
compactness  due  to  the  growth  or  decay  of  ice. 

Ice  strength  is  treated  as  a  function  of  the  ice 
thickness  distribution  and  compactness  given  by  the 
equation 


P  =  P*  h  exp  [  -  C  (1  -  A)  j  , 

where  P*  and  C  are  fixed  empit  ical  constants,  h  is  the 
ice  thickness,  and  A  is  the  compactness.  This  relation¬ 
ship  shows  the  strength  of  ice  to  be  strongly  dependent 
on  the  amount  of  thin  ice  [(l  At).  It  also  allows 
the  ice  to  strengthen  as  it  becomes  thicker. 

The  thermodynamic  portion  of  the  code  determines 
growth  and  decay  rates  of  ice  based  on  a  heat  budget 
balance  between  the  atmosphere,  the  ice,  and  the 
ocean,  including  the  effects  of  heat  absorbed  bv  leads 
via  lateral  mixing.  Similar  to  Semiuer's  (1976) 
formulation,  heat  is  transferred  through  the  ice  by 
assuming  a  linear  temperature  profile  along  with  a 
constant  ice  conductivity.  When  open  water  is  losing 
heat  to  the  atmosphere,  the  heat  budget  growth  rates 
arc  taken  to  be  vertical  growth  rates.  When  open  water 
absorbs  heat,  the  heat  mixes  underneath  the  flows  to 
reduce  the  vertical  growth  rate.  Any  remaining  heat 
can  either  cause  lateral  melting  or  raise  the  tempera¬ 
ture  of  the  mixed  layer.  In  the  presence  of  an  ice  cover. 


the  mixed-layer  temperature  is  always  set  equal  to 
freezing.  Thus,  excess  heat  absorbed  by  leads  is  used 
for  lateral  melting  until  the  ice  disappears.  During 
growth  conditions,  ice  cannot  form  until  the  mixed 
layer  reaches  the  freezing  temperature  of  seawater. 

In  the  two-level  version  of  the  model  (following 
Bryan  et  al.,  1975  and  Manabe  et  ai.,  1979),  the  effects 
of  snow  cover  are  treated  such  that  the  ice  surface 
albedo  is  that  of  snow  (0.75)  when  the  calculated 
surface  temperature  is  below  freezing  and  that  of  snow- 
free  ice  (0.66)  when  the  surface  temperature  is  at  the 
melting  point.  Thus,  the  upward  heat  flow,  Ih, 
through  ice  of  thickness  h  is 

lh  =  (K/h)  <TW  -  TJ  , 

where  K  is  the  ice  conductivity,  Tw  is  the  water 
temperature,  and  Ta  is  the  surface  temperature  of  the 
ice. 

In  the  two-level  case,  snow  is  parameterized  only 
through  the  surface  albedos,  but  the  new  seven-level 
formulism  uses  the  accumulated  rates  from  Maykut 
and  Untersteiner  (1969)  and  Parkinson  and 
Washington  (1979).  The  thermal  conductivitv  in  the 
seven-level  case  is  a  single  value  based  on  a  weighted 
sum  of  snow  and  ice  conductivities 

(K$SnLVL  +  K,hLVL) 

where  hLVL  is  the  ice  thickness  at  that  level,  Sn,  VL  is 
the  snow  depth  at  the  same  level,  A",  is  the  snow 
conductivity,  and  K ,  is  the  ice  conductivity.  The 
prescribed  surface  albedos  used  by  Walsh  et  al.  are  0.80 
for  snow  and  0.65  for  ice. 

The  surface  heat  budget,  after  Parkinson  and 
Washington  (1979)  and  Manabe  et  al.  (1979),  is 
given  by 

(1  -  a)  F\  +  F,  +  D,  |Tg  (Tu  -  TJ 
+  D2  \  U^fqjTj  -  </,  (TJ] 

-  D,  T*  +  (K/H)  (T„  -  TJ  =  0  , 

where  a  is  the  surface  albedo,  T0  is  the  surface  tem¬ 
perature  of  ice,  Ta  is  the  air  temperature,  TK  is  the 
water  temperature,  U  is  the  geostrophic  wind,  qs  is 
the  specific  humidity  of  the  ice  surface,  Fs  is  the 
incoming  short-wave  radiation,  F,  is  the  incoming 
long-wave  radiation,  is  the  bulk  sensible  heat 
transfer  coefficient,  D2  is  the  bulk  latent  heal  transfer 
coefficient  (water  or  ice),  and  D,  is  the  Stephan- 
Boltzman  constant  times  the  surface  emissivity.  This 
surface  iieat  budget  defines  a  surface  temperature  for 
the  ice  that  balances  the  heat  budget.  This  temperature 
then  determines  the  conduction  of  heat  through  the 
ice  and  the  growth  rate.  If  the  derived  temperature  is 
above  freezing,  it  is  set  back  to  the  freezing  point. 


Surface  and  bottom  ablation  rates  are  then  detei  mined 
by  the  imbalances  in  the  surface  heat  budget  and  by 
conduction  of  heat  into  the  mixed  layer.  Heat  transfer 
from  the  deep,  warmer,  ocean  water  can  either  be 
treated  as  a  constant  or  as  a  variable  heat  flux  into 
the  mixed  layer.  For  a  detailed  discussion  of  the 
thermodynamic  portion  of  the  model,  see  Hibler 
(1980). 

III.  Model  Grid 

The  Barents  Sea  model  giid  was  designed  as  a 
subsection  of  the  FNOC  northern  hemisphere  polar 
stereographic  grid.  The  model  grid  covers  the  entire 
Baicnts  Sea  and  the  western  half  of  the  Kaia  Sea.  An 
averaged  mapping  factor  is  used  to  approximate  equal 
spacing  for  the  FNOC  polar  stereographic  grid  in  the 
Barents  Sea  domain.  The  ice  model  grid  is  defined  as 
an  equally  spaced,  25-km  grid,  subset  of  the  FNOC 
northern  hemisphere  polar  stereographic  grid.  The 
resultant  ice  model  dimensions  are  74  ,  66  (Fig.  1). 

Boundaries  of  the  model  are  solid  walls  except  for 
the  boundary  between  Spitzbergen  and  Norway.  This 
region  contains  tw'o  rows  of  “outflow”  grid  cells.  Ice 
can  only  be  transferred  into  these  grid  cells  by  advection 
and  once  there,  flows  out  of  the  basin.  The  boundaries 
between  Spitzbergen  and  Franz  Josef  Land,  Franz 
Josef  Land  and  Novaya  Zemlya,  and  Novaya  Zemlya 
and  the  Soviet  coast  are  all  inflow/outflow  boundaries. 
Details  of  the  model’s  inflow/outflow  boundary 
conditions  will  be  described  in  the  model  development 
section  VII  B. 

IV.  Forcing 

The  Barents  Sea  model  is  driven  by  both  atmospheric 
and  oceanic  forcing.  The  atmosphere  forcing  is 
obtained  from  the  NOGAPS  model.  Ibis  global 
atmospheric  model  provides  surface  pressme  fields 
(NOGAPS  field  A0 1),  which  are  used  to  determine 
geostrophic  winds.  In  addition  to  surface  pressure,  the 
NOGAPS  model  also  provides  surface  vapor  pressure 
(A12),  which  is  used  in  conjunction  with  surface 
pressure  to  determine  the  specific  humidity  at  the  ice 
surface;  surface  air  temperature  (A07);  incoming  solar 
radiation  (short  wave- A  1 1);  sensible  heat  flux  (A  16); 
and  total  heat  flux  (A  18).  These  last  three  fields  are 
used  to  determine  long-wave  radiation. 

Monthly  mean  geostrophic  ocean  currents  and  deep- 
ocean  heat  fluxes  derived  from  the  Hibler/Bryan 
coupled  ice-ocean  model  are  used  as  the  ocean  forcing 
for  the  model.  The  effects  of  the  variability  of  ocean 
currents  on  ice  drift  has  been  shown  to  be  of  importance 
over  long  time  scales  (Thorndike  and  Colony,  1982). 
On  the  time  scale  of  a  forecast  (5  days),  the  variability 
of  the  ocean  currents  has  a  much  smaller  effect  on  the 
ice  drift  than  the  variability  of  the  wind  stress  fields. 
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For  this  reason,  monthly  mean  ocean  currents  can  be 
used  as  a  reasonable  estimate  for  oceanic  forcing. 

Including  monthly  mean,  deep-ocean  heat  fluxes  has 
resulted  in  a  tremendous  improvement  in  the  model’s 
capability  to  predict  edge  location  in  the  marginal  ice 
zone.  Hibler  and  Bryan  (1984;  1987)  have  shown  that 
this  oceanic  heat  flux  can  melt  large  amounts  of  ice 
in  the  marginal  ice  zone.  A  drastic  improvement  in 
ice-edge  location  was  seen  in  the  results  of  the  PIPS 
model  when  monthly  mean,  deep-ocean  heat  fluxes 
were  included  (Preller,  1985;  Preilerand  Posey,  1989). 
The  Barents  Sea  model  required  some  adjustment  to 
the  original  Hibler  and  Bryan  heat  fluxes.  This 
correction  to  the  heat  fluxes  is  discussed  in 
section  VII  A. 

Both  PIPS  and  RPIPS-B  forecasts  could  be 
improved  by  improving  the  atmospheric  and  oceanic- 
forcing .  Oceanic  forcing  would  be  improved  by  hav  ing 
an  ocean  model  coupled  to  the  ice  model.  This  model 
could  predict  the  variability  of  the  ocean  on  the  same 
time  and  space  scales  as  the  ice.  Higher  resolution  in 
the  atmospheric  models  would  also  improve  the 
forecast  capability.  At  present,  NOGAPS  resolution 
is  approximately  400  km,  but  RPIPS-B  is  25  km.  The 
development  of  such  models  is  a  task  being  researched 
in  both  the  (6.1  and  6.2)  communities. 

V.  Initial  Conditions 

During  the  development  of  the  Barents  Sea  model, 
the  model  was  spun  up  to  a  cyclic  “equilibrium  state,” 
where  ice  thickness,  ice  velocity,  and  ice  concentration 
take  on  similar  values  on  corresponding  days  of 
successive  years.  Initialization  for  the  equilibrium  state 
case  requires  setting  ice  drift  velocities  to  zero  and 
ice  thickness  to  a  constant  value  of  1.1  m 
(1.0  x  lO'kgm-Vp,)  and  ice  concentration  to  lUH'u 
at  all  grid  points.  From  these  initial  conditions,  it  takes 
approximately  2  years  of  model  integration  to  reach 
a  cyclic  equilibrium  state.  One  particular  year  of 
NOGAPS  forcing  is  used  repeatedly  to  reach  the 
equilibrium  state. 

RPIPS-B  can  be  initialized  in  three  different  ways. 
Each  day  the  model  is  run  making  a  120-hour  forecast. 
The  model’s  24-hour  forecast  of  ice  thickness,  concen¬ 
tration,  ice  drift,  surface  ice  temperature,  and  heat 
absorbed  by  the  open  ocean  is  saved  each  day.  I  he 
model  uses  this  24-hour  forecast  as  its  restart  field  the 
next  day.  If  the  restart  field  from  the  previous  day  is 
not  available,  then  the  model  searches  back  as  far  as 
I  week.  If  no  restart  fields  can  be  found,  then  model 
climatology  is  used  to  restart  the  model.  The  model 
climatology  contains  monthly  mean  fields  derived  from 
the  Barents  Sea  “cyclic  steady  state”  results  driven  by 
1986  NOGAPS  forcing. 

Once  per  week  (usually  Friday),  in  addition  to  the 
restart  field,  the  model  is  also  given  a  new  laid  of 


gridded  ice  concentration.  This  concentration  field  is 
a  digitized  version  of  the  NPOC  weekly  ice  concen¬ 
tration  analysis.  The  NPOC  analysis  is  a  subjective 
analysis  derived  from  available  remotely  sensed  data 
(Advanced  Very  High  Resolution  Radiometei .  V  isible, 
Passive  Microwave)  and  available  observations  (ship, 
plane,  etc.)  (Naval  Polar  Oceanography  (  enter.  19X6). 
This  field  is  hand-digitized  once  per  week  by  NPOC 
and  transferred  to  FNOC.  The  digitized  data  are  then 
placed  on  the  model  grid  and  used  to  update  the  model. 

The  data  are  assimilated  into  the  Barents  Sea  model 
in  the  following  manner.  The  model’s  restart  field  is 
found  (either  a  24-hour  forecast  oi  climatology). 
RPIPS-B  then  checks  to  see  if  an  NPOC  analysis  field 
is  available.  If  the  analysis  field  is  older  than  4  days 
or  has  been  used  within  the  past  4  days,  then  an  update 
is  not  made.  Otherwise,  if  the  NPOC  analysis  is 
available,  then  the  model's  forecasted  ice  concentration 
field  is  entirely  replaced  by  the  NPOC  analysis.  Two 
additional  fields  are  then  updated:  the  ice  thickness 
and  the  heat  stored  by  the  ocean.  Hie  new  concen¬ 
tration  field  is  compared  to  the  model-derived 
concentration  field.  If  no  ice  exists  where  it  existed 
before  the  update,  then  the  ice  thivkiiess  is  set  equal 
to  zero  and  a  small  amount  of  heat  is  added  to  the 
open-ocean  mixed  layer.  If  conceiitiatiou  lias  been 
added  to  a  previous  open  ocean  icgion,  then  ice 
thickness  is  updated  in  the  following  manner: 

if  0. 15<A  <0.5  and  H<0.2.  then  //  0.2 

or 

If  A  >0.5  and  H<0.2.  then  //  o .4, 

and  heat  is  removed  from  the  mixed  layet. 

VI.  Model  Parameters 

The  Barents  Sea  model  uses  a  6-hour  time-step.  In 
the  initial  model  testing,  atmospheric  forcing  fields 
were  interpolated  from  the  63  x  63  northern 
hemisphere  polar  stereographic  grid  used  by  FNOC  to 
the  Barents  Sea  grid.  Ocean  fotcing  is  updated  once 
per  month.  Turning  angles  for  tC’  g  -ostrophic  wind 
and  drag  coefficients  were  based  on  values  derived 
from  the  testing  of  PIPS  (Preller  and  Posey,  1989). 
Additional  parameters  used  by  the  model  are  given 
in  Table  I . 

VII.  Model  Development 

The  Barents  and  Kara  Seas  are  both  Arctic  maiginal 
seas.  The  Barents  Sea  is  bounded  by  Spiizbergen  and 
Franz  Josef  Land  to  the  north,  and  to  the  west  by  a 
line  running  north  from  North  Cape  tluough  Bear 
Island  to  the  southern  tip  of  Spiizbergen  I  he  Barents 
and  Kara  Seas  are  separated  by  the  island  ol  No'  avw 
Zemlya,  which  restricts  the  transpon  (>!  uc  between 


them.  These  seas  are  both  relatively  shallow;  the 
Barents  Sea,  the  deeper  of  the  two,  averages  200  in  in 
depth. 

The  Kara  Sea  is  typical  for  its  high  northern  latitude. 
It  has  a  solid  ice  cover  during  the  winter  and  is  ice-free 
during  the  summer,  except  for  the  extreme  northern 
part.  As  the  warmest  of  the  Arctic  seas,  the  Barents 
Sea  is  atypical.  It  is  strongly  influenced  by  the  warm, 
saline  Atlantic  waters  of  the  Norwegian  Current.  As 
a  result,  the  southwestern  third  of  the  Barents  Sea 
remains  ice-free  throughout  the  year. 

The  composition  of  the  sea  ice  in  the  Barents  Sea 
is  usually  complex  and  is  composed  of  many  different 
ice  types.  The  majority  of  the  ice  is  first-year  ice  formed 
locally,  but  some  multiyear  ice  is  transported  from  the 
Arctic  Ocean  between  Spitzbergen,  Franz  Josef  Land, 
and  the  northern  tip  of  Novaya  Zenrlya.  Multiyear  ice 
can  also  form  locally  around  some  of  the  islands,  and 
icebergs  frequently  occur  in  the  Barents  Sea.  During 
the  summer,  the  ice  retreats  far  to  the  north  and  the 
Barents  Sea  becomes  ice-free.  In  winter,  the  ice  edge 
advances  until  it  extends  over  the  Spitzbergen  Bank 
to  Bear  Island.  The  ice  often  extends  well  off  the  coast 
of  Novaya  Zemlya  and  over  the  shelf  off  the  Soviet 
coast  (the  Pechora  Sea)  and  into  the  White  Sea. 

Year-to-year  variability  of  the  ice  extent  in  the 
Barents  Sea  is  srrong.  Parkinson  et  al.  (1987;  studied 
Arctic  sea  ice  using  passive  microwave  data  from  the 
Nimbus-5  Electrically  Scanning  Microwave  Radiometer 
(ESMR)  from  1973  to  1976.  They  found  that  the  timing 
of  the  maximum  extent  of  ice  varies  from  year  to  year 


and  is  most  likely  due  to  a  combination  of  atmospheric 
and  oceanic  effects.  Figure  2  shows  their  calculation 
of  the  yearly  cycle  of  sea-ice  covei  troiu  1973  to  1976. 
Note  that  a  maximum  occurs  in  l  ebruan  197  3  and  that 
a  double  maximum  occurs  during  I  ebi  ua. »  and  .\p.  il 
1974.  The  maximum  is  in  March  and  April  1975,  while 
maxima  occur  in  January  and  Maich  1976.  Summer¬ 
time  variability  is  also  apparent  from  tlw  change  in 
the  amount  of  ice  area  existing  in  August  and 
September  from  year  to  year. 

Year-to-year  variability  is  also  seen  in  the  climato¬ 
logical  minimum  and  maximum  extension  of  sea  ice 
in  the  Barents  Sea  from  the  Navy-NOAA  Joint  Ice 
Center  (Fig.  3).  Despite  these  yearls  satiations,  the 
seasonal  pattern  of  growth  and  deca;  in  the  Barents 
Sea  is  similar  from  year  to  sear,  with  the  maximum 
extension  of  the  ice  usually  falling  in  the  Fcbiuary -April 
time  fiame  and  the  minimum  occurring  in  August- 
September  (Parkinson  cl  al.,  1987;  1  oeng,  1979;  l.oeng 
and  Vinje,  1979). 

A.  Ocean  Currents  and  Heal  Fluxes 

The  location  of  the  ice  edge  is  strongls  influenced 
by  the  heat  brought  into  the  Barents  Sea  from  the 
southwest  via  the  Norwegian  Atlantic  Current.  The 
Norwegian  Coastal  Current  Hows  along  the  coast 
of  Norway.  Paralleling  this  coastal  current  is  the 
Norwegian  Atlantic  Current.  Off  the  coao  of  north¬ 
ern  Norway,  this  current  splits  in  two;  one  braneh 
flows  northward  and  becomes  the  West  Spitzbergen 
current;  the  other  branch  enters  the  Barents  Sea  as  the 


Table  1.  Parameters  used  in  RPIPS-B 


Parameter 

Definition 

Value 

ca 

Drag  coefficient  of  air 

0  0008 

Cw 

Drag  coefficient  of  water 

0  U0b5 

c 

Empirical  constant  in  the 
strength  equation 

20 

e 

Ratio  of  the  principal  axis  of  the 
elliptic  yield  equation 

2 

t 

Coriolis  parameter 

1  46  x  10  4  sec 

hO 

ThicSn? limit  between  thick 
and  thin  ice 

0.5  m 

p, 

Density  of  ice 

0  91  x  10J  kg  m  4 

**  a 

Density  of  air 

1  3  kg  m  ! 

p- 

Pressure  constant 

2  75  x  104  N  m  ; 

> 

X 

11 

> 

■< 

Horizontal  grid  spacing 

25  km 

At 

Time  step 

6  hour 

£max 

Nonlinear  bulk  viscosity 

(P/4)  x  109  kg  s  1 

rjmax 

Nonlinear  shear  viscosity 

7ma  x/eJ 

* 

Turning  angle  of  air 

23° 

e 

Turning  angle  of  water 

25° 

H 

Initialization  ice  thickness 

1  1  m 

North  Cape  c  urrent.  I  mure  4  is  a  simplitied  Jcpi.iioti  to  be  a  result  of  the  resolution  ot  ■  I M  l  t,-,  and  lit \ an 

ot  the  surface  current  system  (torn  Midttuu  and  I  ocng  ice  ocean  model.  I  sine  160  km  ei  :.i  i  th.  model 

(1987),  based  on  current  maps  by  l.uusiura  I  I959t  and  cannot  accurately  resolve  smut-  te.it. a s-  >■:  ,  u.  ocean 

Novitskiy  ( 1961 ).  The  southern  part  ot  the  Noah  <  'ape  <  iirrent  field.  As  a  result,  the  .  uncut  tic!,!  ■>  smoothed 
Current  continues  ea  .aid  with  the  Norwegian  and  can  be  weaker  than  ohset \  ed.  Weak  ,  u:  on s  would 

Coastal  Current  arm  ‘  ecotnes  the  Mur  man  Cut  t  cut .  not  penetrate  far  enough  i  ast  ,n;,  tic  I;  ,m:t.  Sea  to 

The  northern  pc  of  the  North  Cape  (intent  splits  bring  the  heat  from  the  Not  wee. a. t  .'.irio,,  ne  cssary 

into  three  brunches.  One  branch  proceeds  northwaid  to  keep  the  region  ice-tree.  I>>  son  \t  'ot  t!. i -  co  essive 
between  Hopen  Island  and  the  Great  Bank  until  it  amount  of  ice.  the  heat  tluxcs  wcu  .  i  a  to  .  o' er 
subr.ierges  under  Atlantic  water.  The  second  branch  a  larger  part  of  the  eastern  harem-  '*  I  h.  >c  heat 


■  mtinues  east w aid  between  the  Cuc.it  Bank  and  the  fluxes  were  adjusted  so  that  ti  e  me  i.  .  ,  e-dee  location 
Central  Bank  as  an  intermediate  e intent  I  Its  thud  agreed  closely  wait  the  NT»  h  ",u  ,r  :ce  edge, 
branch  turns  southeastward  and  parallels  the  Mia  nan  New  heat  fluxes  were  created  f-v  -a  ,_n.g  idtacent 

current  oft  the  coast  ot  the  Soviet  Union  and  tits  n  tiiitts  row  (column)  values  and  uisoii,;--  ti.  :t- ■■■  oerage 

northeastward  along  the  axis  of  the  eastern  basin  oft  values  between  the  old  row s  icolunm- 1  j !  w  ,i>  done 

the  coast  ot  Nov aya  /emlya.  The  colder,  less  -.aline  in  such  a  way  that  the  large  mud  l„i\  aim-  were 

Arctic  water  enters  the  sea  as  two  currents,  one  along  expanded  and  the  low  values  close  to  land  boundaries 


the  east  coast  ot  Spit/bergen  and  one  between  1  tan/  were  eliminated.  Figures  8a  and  Sh  show  the  improve- 
Joset  Land  atid  Nov  aya  /.emlya  (Dickson  et  ah,  1970).  ment  to  the  ice  edge  location  when  the  model  was 

The  first  current,  the  hast  Spit/beigen  current.  driven  by  1983  NOGAPS  forcing  and  the  adjusted  heal 

continues  southward  along  the  coast  of  Spit/bergen.  fluxes.  Figures  9a  and  9b  show  the  imptov ement  in 

The  second  current ,  the  Persey  C  ut  rent,  flows  south-  model  results  when  the  model  was  (.haven  by  1986 

westward,  south  of!  ran/ Josef  land,  along  tiie  eastern  NOGAPS  forcing  and  the  adjusted  heat  tluxes.  The 

slope  of  the  Svalbarvl  Bank  and  around  Bear  Island.  correction  to  the  heat  fluxes  was  chosen  as  the  best 

Examples  ot  the  ocean  currents  used  in  the  Barents  average  adjustment  between  results  using  the  two  years 

Sea  model  are  shown  in  figures  5a  and  5b.  As  stated,  of  atmospheric  forcing, 

these  monthly  mean  values  are  derived  from  the  Ilibler 

and  Bryan  ice-ocean  model.  1  hese  currents  follow  the  B.  Open  Boundary  C  onditions 
general  pattern  described  above.  However,  due  to  In  the  initial  testing  of  the  Barents  sea  model,  it 

the  coarse  resolution  of  the  Htbler  and  Bryan  model  became  obvious  that  the  region  required  ice  inflow 

( =  1 60  km),  some  ot  the  detail  of  the  cut  rent  pattern  out  How  boundary  conditions.  Prior  to  this,  both  PIPS 

is  lost,  |n  most  of  the  monthly  mean  ctments.  the  and  tests  performed  by  Hihler  had  med  only  ice 

inflowing  Atlantic  water  breaks  up  into  two  as  opposed  outflow  boundary  conditions. 

to  thiec  branches;  one  moves  to  the  east  along  the  The  first  test  of  the  model  hoiiiidar.  treated  all 

Norwegian  coast .  and  a  second  enters  the  B.uents  Sea  boundaries  as  solid  walls  except  tot  the  boundary 

and  turns  toward  the  not  (beast.  Dtirutg  must  mouths,  between  Spit/bergen  and  Norway .  1  hat  bound. try  is 

the  current  enters  the  Barents  Sea  between  I  i.m/  Josef  treated  as  an  ice  outflow  boundary  in  all  ot  the  te't 

1  and  and  Nov  ay  a  /emlya.  flows  low  at  J  the  southwest.  cases,  lee  can  be  I  rails  fared  into  these  id  a-li'  only 

and  curves  hack  tow.ud  the  north.  Par.  >-l  t ! - : - ,  intent  by  advection  and,  once  theic.  is  umow-i  i-ou,  the 

exits  through  t tic  teeiott  botnet  n  i  tan/  f  ,  -  i.i.md  hasin.  Allowing  outflow  only  ti  (h.s  i !ai  i- .»  valid 

Spit/bergen,  and  pai!  I -vs  aloi.-c  Spit  -I  u  .  n-l  «-tt  assumption,  since  most  of  the  How  ts  tnt*»  the  Barents 

ot  the  Barents  Sea  nwr  Bear  Island  H  tg.  'a;  I  lining  Sea  on  the  southern,  ice -free  half  of  the  houndaiy  and 

October  and  November,  howeva .  the  lie’*  at-  the  out  of  the  Barents  Sea  on  the  iioithciii  -ule  ol  the 

Barents  Sea  through  the  passage  bet w ecu  I  t .-.it/  loset  boundary  where  ice  can  exist ,  It  should  .d  o  tie  noted 

l  and  and  Nov  ay  a  /emlya  (fig.  5b).  that  these  early  tests  were  pa  lot  med  o.,  a  version  of 

The  heat  fluxes  associated  with  these- cm  tents  tot  the  model  that  excluded  the  Kata  Ss.i.  When  the 
the  Htbler  and  Bivan  modeled  region  ate  shown  in  northern  boundaries  of  the  model  vute  kept  closed. 

Figure  6.  I  his  result  indicates  that  l.ng.  In  Huxe-  the  model  gave  very  unrealistic  icsiih  It  the  wind 

are  carried  into  the  mu.Hi  sc-;.-tn  »ln:d  -t  . I  ■  i.  •  „  forcing  blew  act oss  these  bound. m.  -  i , .  m  "te  north. 

Sea.  The  mom  hi  v  mean  heat  (luxes  !  n .  :  o  >i:e  then  ice  was  pushed  away  and  lug.  ii.nl  ot  regions 

Barents  Sea  grid  ate  shown  in  I  iguiis  ,,  of  open  water  appeared  at  iltc  solid  Toutidarv . 

f  igure  7a  is  the  original  monthly  mean  heat  ibo  I  eld  If  the  wind  forcing  was  trom  the  somli.  dien  ice  would 

from  the  ice-ocean  model,  \fter  cxfe.isi*  c  u-sm.g  of  “pile  up"  unrealistically  at  the  bounda.v  (l  ie  10a). 

the  model  w  ith  these  heat  tluxes,  n  was  detcniimc.l  that  In  the  second  ease,  the  hound. h  v  Ktc-ccd  Iran/ 

the  (luxes  were  too  weak  m  the  south  ennui  ution  Josef  1  and  and  Novaya  /eml  ,a  vvm  ,k  -ilmi  md  an  ice 

ot  the  Barents  Sea.  \s  a  result,  tin-  teg  nut  nlv.avs  “inflow  outflow"  boundat  . .  I  Ik  ti  • ...  v  outflow 

contained  too  much  t  e  I  Itese  weak  heat  tluxes  ap|  eat  boundary  is  the  first  intetiot  ".oumat  -t  gt  id  cells 


between  Franz  Josef  Land  and  Novaya  Zemlya.  For 
each  time  step,  the  model  checks  the  ice  drift  at  each 
grid  cell  on  the  boundary.  If  drift  implies  that  ice  is 
flowing  out  of  the  basin,  then  the  standard  Hibler  ice 
outflow  boundary  'ondition  is  used.  If  the  drift 
indicates  inflow,  a  specified  amount  of  ice  thickness 
is  placed  into  the  grid  cell  and  is  allowed  to  advect  into 
the  model.  The  amount  of  ice  placed  into  the  grid  cell 
is  determined  in  one  of  two  ways.  In  the  earlv  testing 
of  the  model,  the  value  of  inflowing  ice  was  set  to  a 
constant  0.5  m  for  all  months  of  the  year,  except 
August  and  September,  when  it  is  set  equal  to  0.0  m. 
During  these  2  months,  the  region  north  of  the  Barents 
Sea  is  most  often  ice-free.  This  value  proved  to  be  a 
good  average  for  ice  inflowing  at  that  boundary.  In 
RPIPS-B  and  additional  model  test  cases,  the  ice 
thickness  derived  from  the  PIPS  model  is  interpolated 
to  the  Barents  Sea  model  grid  and  used  as  the  ice  inflow 
value.  Although  the  second  method  is  more  accurate, 
the  constant  value  used  in  the  first  method  was  very 
useful  m  that  it  allowed  the  Barents  Sea  model  to  be 
run  as  a  "stand-alone”  model.  That  is,  it  allowed  many 
test  cases  to  be  run  when  the  PIPS-intcrpoiated  ice 
thickness  values  were  not  available.  Figure  10b  shows 
a  result  identical  to  that  of  Figure  10a,  except  that  the 
boundary  between  Franz  Josef  Land  and  Novaya 
Zemlya  has  been  opened  and  a  0.5-m  average  ice  thick¬ 
ness  applied  at  the  inflow/outflow  boundary.  This 
result  was  in  much  better  agreement  with  the  NPOC 
analysis  than  the  closed  boundary  result. 

In  the  third  test  case,  the  same  inflow /outflow 
boundary  condition  was  added  to  the  boundary 
between  Franz  Josef  Land  and  Spitzbergen.  Figure  11a 
shows  the  result  when  this  boundary  was  closed  and 
Figure  1  lb  is  the  identical  case  result  with  the  boundary 
open  to  inflow  and  outflow.  Again,  opening  this 
boundary  resulted  in  a  much  more  realistic  solution. 

In  the  RPIPS-B  model  and  the  remaining  test  cases, 
the  western  Kara  Sea  was  added  to  the  model.  In  the 
Kara  Sea,  however,  the  constant  value  of  0  5  m  is  used 
as  the  boundary  condition.  This  method  was  chosen 
because  the  PIPS  grid  does  not  cover  all  of  the 
Kara  Sea. 

The  remaining  test  cases  discussed  in  this  ieport  and 
the  RPIPS-B  model  use  open  inflow/outflow  bound¬ 
ary  conditions  from  PIPS  at  the  boundaries  between 
Spitzbergen  and  Franz  Josef  Land  and  between  Franz 
Josef  Land  and  Novaya  Zemlya,  and  a  constant  value 
inflow/outflow  boundary  condition  between  Novaya 
Zemlya  and  the  Soviet  coast  (Kara  Sea).  The  western 
boundary  between  Spitzbergen  and  Norway  is  treated 
as  an  ice  outflow  boundary  condition  in  all  cases. 

VIII.  Example  Output  from  RPIPS-B 

Results  from  the  Barents  Sea  model  test  cases  using 
1983  and  1986  NOGAPS  atmospheric  forcing,  along 


with  results  from  the  RPIPS-B  operational  check¬ 
out  (OPCHECK),  will  be  presented  in  this  section. 
During  OPCHECK,  the  model  johstrcam  is  tun  to 
show  that  the  job  can  perform  m  an  operational 
environment  without  disruption  of  the  operational  mu. 
It  also  provides  an  opportunity  for  minor 
improvements  to  be  made  to  the  model.  The  main 
difference  between  the  Barents  Sea  model  test  results 
and  the  RPIPS-B  results  is  that  RPIPS-B  results  are 
updated  weekly  by  the  NPOC  analysis  of  ice  concen¬ 
tration.  Model  results  will  be  examined  in  terms  of  ice 
drift,  ice  thickness  and  ice  concentration  (ice  edge). 

A  number  of  similar  tendencies  appear  in  Barents 
Sea  model  solutions  using  both  the  1983  and  1986 
NOGAPS  forcing.  Therefore,  we  will  present  results 
from  the  evolution  of  the  1983  test  case  to  imply  that 
a  similar  situation  existed  in  the  1986  case.  The  Barents 
Sea  model  was  shown  to  take,  at  most,  two  years  to 
spin  up  to  a  cyclic  equilibrium  solution.  Figures  12a 
and  12b  show  the  March  mean  thickness  field  for  the 
first  and  second  year  of  spin-up.  The  model  1983  test 
run  was  initialized  on  January  18  using  forcing  from 
NOGAPS  and  the  initial  conditions  described  in 
section  V.  The  model  is  integrated  out  for  3  years  using 
1983  NOGAPS  forcing  to  drive  the  model  each  year. 
The  main  difference  between  the  first-year  and  second- 
year  results  appear  in  the  ice  thickness  fields  near 
Spitzbergen  and  Franz  Josef  Land.  These  tesuhs  imply 
that  the  model  needs  longer  than  2  months  of  spin-up 
time  (Fig.  12a)  to  build  up  thick  ice  along  these  islands 
(Fig.  12b).  The  third  year  monthly  mean  ice-thickness 
field  is  identical  to  the  second-year  ice  thickness  and 
is  not  shown. 

Figure  13  shows  the  RPIPS-  B  grid  with  four  transect 
lines  overlaid  at  points  J  --  5,  20,  35  and  50. 
Figures  14a-d  show  March  monthly  mean  ice  thickness 
along  each  transect  for  the  first  and  second  sear  of 
model  integration.  The  largest  difference  from  year  1  to 
year  2  lies  along  the  transect  at  J  --  5,  Figure  14a.  This 
transect  runs  along  the  eastern  side  of  Spitzbergen  and 
shows  the  piling  up  of  ice  along  this  coast  during  the 
second  year  winter.  Transects  from  the  third  year  of 
integration  are  identical  to  the  second-year  transects 
and  are  not  shown. 

The  Barents  Sea  is  observed  to  base  its  most 
extensive  ice  cover  in  the  w  inter  months  (February  and 
April)  and  to  become  ice-free  or  nearly  ice-free  in  the 
summer  (August  and  September).  Model  results 
develop  a  similar  seasonal  pattern,  w ith  the  Barents  Sea 
becoming  ice-free  during  the  first  summer  of  the 
spin-up  testing.  Results  show  that  after  the  first  year's 
summer  melting  of  ice  the  model  quickly  reached 
a  near-equilibrium  state  by  fall  of  the  first  year. 
Figures  15a  and  15b  are  November  mean  ice  thicknesses 
from  the  first  and  second  year  of  mode!  integration. 
Only  small  differences,  again,  in  the  region  between 
Spitzbergen  and  Franz  Josef  1  and,  can  be  seen. 


Figures  16a-d  show  November  mean  thicknesses  along 
the  four  transects  indicated  in  Figure  13.  Aside  from 
the  small  increase  in  ice  thickness  at  transect  .1  -  5, 
these  figures  show  little  change  in  ice  thickness  from 
year  1  to  year  2. 

Equilibrium  results  from  both  the  1983  and  1986  test 
cases  showed  that  the  thickest  ice  in  the  Barents  Sea 
occurs  in  the  winter  months  on  the  eastern  coast  of 
Spitzbergen,  where  model  thicknesses  are 
approximately  4  m.  Throughout  the  rest  of  the  Barents 
Sea,  the  ice  thickness  averages  about  1  m.  During  the 
summer  months,  when  ice  is  present,  it  averages  less 
than  1  m  thick. 

Ice  drift  in  the  Barents  Sea  is  closely  correlated  to 
the  winds.  The  ice  drift  from  the  Barents  Sea  model 
responds  rapidly  to  the  wind  forcing.  This  rapid 
response  to  the  winds  has  been  observed  in  the  Barents 
Sea  (Vinje,  1985). 

Ice-thickness  values  in  many  parts  of  the  Barents  Sea 
appear  to  be  closely  linked  to  the  wind  forcing. 
Figures  17a  and  17b  show  April  monthly  mean  ice  drift 
and  ice  thickness  from  1986.  Monthly  mean  winds 
during  this  period  were  blowing  toward  the  south  and 
the  southeast.  As  a  result,  ice  piled  up  on  the  western 
side  of  Novaya  Zemlya  and  along  the  Soviet  coast  in 
the  Kara  Sea.  Figures  18a  and  18b  show  that  in  May, 
the  mean  winds  changed,  causing  the  ice  to  drift  away 
from  the  Soviet  coast,  awav  from  the  west  coast  of 
Novaya  Zemlya,  and  into  the  east  coast  of  Spitzbergen. 
Correspondingly,  the  ice  thinned  along  the  Soviet  coast 
and  the  west  coast  of  Novaya  Zemlya  and  piled  up 
along  the  east  Spitzbergen  coast. 

The  high  resolution  of  the  Barents  Sea  model  allows 
for  realistic  simulations  of  ice  moving  toward  or  away 
from  model  coastlines  that  are  not  possible  in  the 
lower-resolution  PIPS  model.  Regions  of  thin  ice  or 
open  water,  induced  by  wind  forcing,  are  often 
observed  in  satellite  imagery  (Fig.  19)  taken  along  the 
coastlines  in  the  Kara  Sea  and  along  the  coasts  of 
Novaya  Zemlya.  The  Barents  Sea  model  is  capable  of 
reproducing  similar  wind-induced  features  in  the  ice 
thickness  field. 

Monthly  mean  ice  concentrations,  the  percentage  of 
a  grid  cell  covered  by  thick  ice,  are  shown  in  Figmes  20a 
and  20b  for  l°8f  and  in  figures  21a  and  21b  for  1986. 
The  NPOC  ■  Jge  is  overlaid  on  the  winter  and  fall 
example:.  1  .arents  Sea  model  tests  show  that  the 
•^odei  res^..s  generally  agree  with  observations  at  all 
times  of  yeai  e  e  t  during  the  summer.  The  Barents 
Sea  moHel  bee. •••  es  ice-free  too  early  in  the  summer 
(usually  miu-  to  late  July)  and  grows  ice  back  too  late 
in  the  fall  (mid-  to  late  October).  The  NPOC  analysis 
shows  the  Barents  Sea  becoming  ice-free  from  mid-  to 
late  August  and  shows  ice  growing  back  in  early 
October.  This  indicates  that  the  Barents  Sea  model 
becomes  too  warm  in  the  summer.  The  extra  heat  that 


keeps  the  Barents  Sea  model  ice-tree  comes  from  a 
combination  of  heating  by  both  the  atmosphere  and 
the  ocean.  This  problem  could  K  «.oi  levied  Mth 
improved  atmospheric  forcing  in  the  Aictic  and 
improved  ocean  forcing  via  a  coupled  ice-ocean  model. 
These  improvements  are  presently  being  addtessed  by 
FNOC  and  NORDA,  and  the  problem  is  being 
corrected  RPIPS-B  by  updating  the  model  with  the 
NPOC  analysis.  This  update  places  ice  in  the  model, 
according  to  the  NPOC  analysis,  and  removes  the 
excess  heat. 

Year-to-year  variability  in  sea-ice  extent,  similar  to 
that  seen  by  Parkinson  et  al .( 1 987 )  is  apparent  in  the 
model  results.  Figure  22  is  a  calculation  of  the  yearly 
sea-ice  cycle  from  the  Barents  Sea  test  cases  using  1983 
and  1986  NOGAPS  forcing.  In  1983,  a  double 
maximum  in  ice  extent  occurred  in  February  and  in 
April  and  May,  similar  to  the  Parkinson  et  al.  results 
for  1974  (Fig.  2b).  The  1986  Barents  Sea  model  results 
showed  a  peak  in  March  and  a  decrease  in  April  and 
May,  similar  to  the  1976  ESMR  data  (Fig.  2d).  Model 
results  also  indicated  a  greater  ice  extent  in  1986  than 
in  1983  (1.4  x  106km-  vs.  1.2  x  10hkm:).  NPOC 
analyses  for  these  2  years  support  this  trend  except  in 
April  1986.  The  NPOC  analysis  showed  ice  cover  as 
extensive  in  April  as  it  was  in  March.  These  figures 
also  demonstrate  that  the  Barents  Sea  model  melts  too 
much  ice  in  the  summer  of  both  years.  Ice  melted  earlier 
in  July  1983  than  in  July  1986,  but  also  grew  back  more 
quickly  in  November  1983.  This  trend  also  appeared 
in  the  NPOC  analysis. 

RPIPS-B  began  to  run  in  an  OPCHECK  mode  in 
the  fall  of  1988.  During  the  model  OPCHECK,  the 
following  fields  were  used  for  testing: 

1.  Cumulative  ice  drift — Tau  24,  Fan  120 

2.  Ice  thickness — Tau  24,  Tau  120 

3.  Ice  concentration — Tau  24,  Tau  120. 

It  is  planned  that  the  operational  version  of  RPIPS-B 
will  generate  the  same  product  fields  as  PIPS  (Preller 
and  Posey,  1989). 

As  stated  previously,  the  only  difference  between 
RPIPS-B  and  the  Barents  Sea  test  eases  is  that  the 
model  is  updated  weekly  with  ice  concent  i  at  ion  data 
from  NPOC.  Both  PIPS  and  RPIPS-B  are  updated 
by  the  same  NPOC  analysis  each  week.  The  finer 
resolution  of  RPIPS-B,  however,  allows  for  a  more 
detailed,  accurate  update  of  the  ice  edge  in  the  Barents 
Sea.  Figure  23  is  the  NPOC  analysis  for  28  December 
1988.  Figures  24a  and  24b  show  the  model  solutions 
after  an  update  from  this  NPOC  analysis  for  both  PIPS 
and  RPIPS-B.  Note  that  the  higher  resolution  and 
larger  areal  extent  of  RPIPS-B  is  responsible  for  the 
detailed  ice  edge  along  Novaya  Zemlya  and  the  Sov  iet 
coast.  The  resolution  of  RPIPS- B  also  allow  s  foi  better 
definition  of  the  ice  protrusion  located  neat  73  N  and 
40°  E 
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IX.  Summary  and  Conclusions 

RPIPS-B,  the  FNOC  version  of  the  Barents  Sea 
model,  is  presently  in  its  OPCHECK  phase  being 
prepared  for  a  winter/spring  OPTEST.  In  its  develop¬ 
ment  stage,  the  Barents  Sea  model  was  spun  up  to  a 
cyclic  equilibrium  state  using  atmospheric  forcing  from 
NOGAPS  and  oceanic  forcing  from  the  Hibler  and 
Bryan  ice-ocean  model.  New  ice  inflow'/outflow 
boundary  conditions  were  designed  for  the  northern 
boundaries  of  the  model. 

Model  test  results  compared  well  with  observations 
(NPOC  analysis).  Results  also  showed  that  the  finer 
resolution  of  RPIPS-B  enabled  far  better  predictions 
than  PIPS  of  ice-edge  location,  ice  thickness,  and  ice 
movement  in  the  Barents  Sea,  particularly  near  land 
boundaries. 

As  a  forecast  tool,  RPIPS-B  will  be  run  daily, 
making  a  120-hour  forecast.  The  model  will  be  run 
after  PIPS  and  will  use  the  ice  thickness  field  from 
PIPS  at  its  inflow/outflow  boundary.  RPIPS-B  will 
produce  the  same  14  product  fields  presently  produced 
by  PIPS. 
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Figure  1.  Barents  Sea  model  grid  with  arrow  indicating  outflow  and  inflow /outflow  boundaries.  Inflow /outflow 
and  outflow  grid  cells  are  indicated  by  "X.  ” 


1 1 


PSEUDO  SEA  ICE  AREA  (106  Km2)  PSEUDO  SEA  ICE  AREA  (106  km2 


Figure  2.  Yearly  cycle  for  the  Barents  Sea  area  covered  by  sea  ice  for  1973-1976  (from  Parkinson  et  al. .  1987).  The 
curves  are  the  ocean  area  covered  with  ice  greater  than  15 %  (15%),  35%  (P35),  50%  ( P50 ),  65%  ( P65 ),  85%  (P85). 
The  pseudo-actual  ice  area  is  the  integrated  ocean  surface  area  covered  by  sea  ice,  excluding  leads  and  polvn  vas. 
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Figure  3.  Climatological  maximum  ami  minimum  ice-edge  locations  from  the  Navy /NO  A  A  Joint  Ice  Center 
for  the  months  of  (a)  March,  (b)  June,  (cj  September  and  (d)  December.  The  dashed  line  is  the  maximum  and 
the  solid  line  the  minimum. 


Figure  4.  Surface  currents  system  for  the  Barents  Sea  from  Midttun  and  L.oeng  (1987). 
Solid  arrows  are  A  tlantie  currents,  dashed  arrows  are  A  retie  currents,  and  dotted  arrows 
are  coastal  currents. 


Figure  5.  Monthly  mean  ocean  current  forcing  for  the  Barents 
Sea  model  derived  from  the  Hibler  and  Bryan  (1987)  ice-ocean 
model.  Monthly  means  for  (a)  February  and  (b)  November 
are  shown.  Maximum  vector  is  5  cm  /sec.  Every  fourth  vector 
is  plotted. 
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Figure  6.  Average  annual  heat  gained  by  the  upper  layer  of  the  ocean  from  the  deep,. ,  <>,  can  ami  lateral 
transport.  Contours  are  in  capacity  of  heat-to-melt  meters  of  ice  per  year  (l  m  year  9.57  U  '-m:) 
(from  Hibler  and  Bryan,  1987). 
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Figure  7.  February  monthly  mean  heat  fluxes  in  watts 
per  meter -  interpolated  to  the  Barents  Sea  grid  from 
(a)  the  Hibler  und  Bryan  model  and  (b)  adjusted  Hibler 
and  Bryan  values. 
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figure  10.  JuU  I  983  ne  thickness  fields  for  (a)  solid  wall 
boundaries,  except  for  an  outflow  boundary  between 
Spitzbergen  and  \orway ;  and  (b)  the  same  case,  except 
that  the  boundary  between  /  fanz  Josef  l. and  and  Novaya 
Zemlya  is  open  to  inflow  and  outflow.  Contour  interval 
is  0.5  m. 


Figure  13.  The  Barents  Sea  grid  with  four  transects  overlaid  at  J  =  5,  20,  35,  50. 
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Figure  15.  November  1983  monthly  mean  ice  thicknesses 
from  (a)  the  first  year  of  model  spin-up  and  ( b )  the  second 
year.  Contour  interval  is  0.1  m. 
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Figure  17.  April  1986  monthly  mean  for  (a)  ice  drift  and 
(b)  ice  thicknesses.  Maximum  vector  for  the  ice  drift  is 
50  cm /sec.  Every  fourth  vector  is  plotted.  Contour  interval 
for  ice  thickness  is  0.5  m. 
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ice  thicknesses.  Maximum  vector  for  the  ice  drift  is 
50  cm/sec.  Every  fourth  vector  is  plotted.  Contour  interval 
for  ice  thickness  is  0.5  m. 


Figure  !9.  NOAA-9  A  VHRR  channel  4  imagery  for  the  Kara  Sea  on  21  February  1986.  The  eastern  side 
of  Novaya  Zemlya  has  new  thin  ice,  indicating  that  ice  was  recently  pushed  away  from  this  coastline  (figure 
courtesy  of  Jeffrey  Hawkins,  NORDA). 
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Figure  21.  The  1986  ice-concentration  fields  for  (a)  March 
and  (b)  November.  Contour  interval  is  in  tenths  or  10%. 
The  ice  edge  from  the  Navy/NOAA  Joint  Ice  Center  is 
overlaid. 
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igure  23.  NPOC  analysis  for  the  eastern  Arctic  on  28  December  1988. 


f  igure  24.  Concentration  field  from  29  December  1988, 
immediately  after  an  update  from  (a)  PIPS  and  (h)  RPIPS-B. 
Contour  interval  is  (a)  5%  and  (h)  10%. 
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